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\Q ' have a de Sitter-invariant ground state for a minimally-coupled massless scalar 
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Tj- , field that gives definite expectation values for any observables not invariant 

under constant shifts of the field. However, if one restricts to observables that 

are shift invariant, as the action is, then there is a unique vacuum state. Here 

j> ■ we calculate the shift-invariant four-point function that is the vacuum expec- 
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^ , tation value of the product of the difference of the field values at one pair of 
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points and of the difference of the field values at a second pair of points. We 
show that this vacuum expectation value obeys a cluster-decomposition prop- 
erty of vanishing in the limit that the one pair of points is moved arbitrarily 
far from the other pair. We also calculate the shift-invariant correlation of 
the gradient of the scalar field at two different points and show that it also 
obeys a cluster-decomposition property. Possible relevance to a putative de 
Sitter-invariant quantum state for gravity is discussed. 



* Alberta-Thy-6-12, larXiv: 1204.44621 [hep-th] 
^Internet address: profdonpage@gmail.com 
■^Internet address: xwu5@ualberta.ca 



1 



1 Introduction 



A quantum state may be denned by the expectation values it gives for quantum 
operators. It is well known that for calculating the expectation values of all products 
of fields at distinct points and of suitably smeared field operators, there is no de 
Sitter-invariant vacuum state for a massless minimally coupled scalar field in a de 
Sitter background [H El El IH EJ E] . Indeed, taking the massless limit of the two- 
point function in the de Sitter-invariant Bunch-Davies (or Euclidean) vacuum of a 
free massive scalar field leads to a divergence, so that this expectation value cannot 
be both a well-defined finite value and de Sitter invariant. 

In order to obtain an infrared (IR)-finite result for the two-point function of a 
massless scalar field, one generally has to abandon the de Sitter invariance of the 
vacuum. For example, the two-point function may have time-dependent terms which 
break de Sitter invariance [5j [6] . 

Since the action of the massless minimally coupled scalar field involves only 
derivative terms, it has a global symmetry under an arbitrary constant shift in field 
values, cf) — > cf) + const. We propose that one should impose this symmetry as a 
physical requirement on the operators whose expectation values define the quantum 
state, thus excluding the non-shift-invariant operators whose expectation values are 
infinite or undefined in the massless limit of the de Sitter-invariant Bunch-Davies 
vacuum. In this paper we shall call non-shift-invariant operators unphysical, so 
physical observables shall here be required to be shift invariant. 

For example, physical n-point functions should be restricted to those that can 
be written in terms of differences of fields at different points and/or in terms of field 
derivatives. Operators like (j)(x)<p(y), {<fi(x),(f)(y)} or (fi 2 (x) — (f) 2 (y), which are not 
shift invariant, do not correspond to physical observables for a massless scalar field 
with a shift-invariant Lagrangian, but operators like [4>{u) — 4>(v )][(f)(x) — 4>(y)] do. 
Then, as we shall see, the above-mentioned IR divergence is only associated with the 
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unphysical correlation functions. For correlation functions made of shift-invariant 
operators, the IR divergences cancel out. A similar viewpoint has been expressed 
by Kirsten and Garriga [7j, who obtained several results overlapping ours. Related 
constructions for the massless scalar field have also been carried out using Krein 
spaces by Bertola, Corbetta, and Moschella [8] and by Bros, Epstein, and Moschella 

m- 

Diffeomorphism invariance for the graviton h a b takes the form of the gauge in- 
variance of its action under 

hab ->• hab + V a £& + V 6 £ a , (1) 

where £ a is an arbitrary infinitesimal coordinate transformation. The graviton two- 
point function is IR divergent in some physical (e.g., transverse-traceless and syn- 
chronous) gauges in the spatially flat patch of de Sitter [10]. However, it is not gauge 
invariant in the sense of (PQ). Indeed, it was shown that one can gauge away the IR 
divergence in the physical gauge [TT] and in some covariant gauges [12] . Moreover, it 
is shown in [T3] that in the open de Sitter patch the two-point function in the phys- 
ical gauge is free of such IR problems. These results suggest that the IR divergence 
of the graviton propagator may not contribute to the physical two-point functions, 
i.e. those constructed from gauge-invariant operators, such as the product of two 
linearized Weyl tensors [HJ [T5J [H] . 

The massless scalar field with the shift invariance considered here is an analogue 
of the graviton with gauge invariance. Indeed they share the same equation of motion 
when the graviton is in the physical gauge [17] . As we shall see, the IR divergence of 
the massless scalar field can be removed if the shift invariance is taken into account, 
so that one only requires the quantum state to give expectation values to physical 
shift-invariant operators. Then correlation functions of physical operators are IR 
finite. We hope to shed some light on the graviton in de Sitter by investigating its 
massless scalar field counterpart as a poor-man's model for the graviton. 
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2 Useful Properties of the de Sitter Geometry 
and Scalar Fields 



Here we consider 4- dimensional de Sitter spacetime. The generalization to other 
dimensions is straightforward. We use units in which the asymptotic Hubble con- 
stant H = a/A/3 is unity for simplicity of notation. In other units, one may use 
dimensional analysis to restore the correct powers of H. Then 4D de Sitter can be 
described by the hyperboloid in 5D Minkowski space, 

VAB X A (x)X B (x) = 1, (2) 

where we use X to denote Minkowski coordinates, and we use x for points on de 
Sitter. More details about the following properties can be found in [18j [19] . 

The geodesic distance /i in de Sitter is, with the affine parameter A, and an 
overdot representing a derivative with respect to that affine parameter, 

fx(x,x') = J d\\/ g a bX a x b . 

A more convenient quantity equivalently characterizing the geodesic distance is de- 
fined by 

Z(x,x') =t 1ab X a (x)X b (x') } (3) 

which is related to \x by 

fj,(x, x') = cos" 1 Z, (— 1 < Z < 1), 

jj(x, x 1 ) = cos _1 (Z — ie), (Z > 1), (4) 

where /i is real for spacelike geodesies (— 1 < Z < 1) and imaginary for timelike ones 
(Z > 1), where fi = ir in terms of the geodesic proper time separation r. Note that 
there is a branch cut along the real axis for Z > 1 where the values of /x on both 
sides are pure imaginary and differ by a sign. Here we choose the convention to pick 
out the limit from below the real axis. Z < — 1 corresponds to the case where the 
separation between the two points is spacelike, but there is no geodesic connecting 
them. In this case \x can still be defined via analytic continuation [T8l [T9] . 
There are two useful unit vectors defined by 

n a = V a /i(x, x'), n v = S/ b 'ii(x, x'), 
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which are pointing outward at each end of the geodesic. They should be distin- 
guished from the normally defined tangent vectors of the geodesic, which are always 
pointing in the direction of increasing A. In particular, for timelike geodesies, n is 
imaginary. So n a n a = n b ' riy = 1 for both spacelike and timelike cases. Note also 
that 

gab>n b ' = -n a , (5) 

V a n 6 > = -(esc ii)(g a y + n a n b >), (6) 

where g a y is the parallel propagator which parallel transports a vector along the 
geodesic. 

In this paper we focus on the Wightman two-point function, defined by the 
vacuum expectation value 

G w (x,y)^(0\<j ) (x)<j ) (y)\0} J (7) 

for a minimally-coupled scalar field of mass m, which obeys the homogeneous equa- 
tion of motion 

(D x -m 2 )G w (x,y) = 0. (8) 

Various Green's functions, such as the Hadamard function or Feynman function, 
can be obtained from G w (see, e.g. [20]). 

For de Sitter-invariant states, this two-point function can only depend on the de 
Sitter-invariant distance, so G w (x,y) = G w (Z(x,y)). If one introduces 

z(x,x') ee l ~ Z ^ x ') = 1 - r]AB [X A (x)-X A (x')}[X B (x)-X B (x')} = sin 2 \^x% 

(9) 

which is one-fourth the square of the distance between the points in the 5D Minkowski 
spacetime and for nearby points is also approximately one-fourth the square of the 
geodesic distance fi(x,x') in the 4D de Sitter spacetime itself, then the Wightman 
two-point function obeys the equation 

z(l - z)f-G w + (2 - 4z)4-G w - m 2 G w = 0. (10) 
dz 2 dz 

In general, the de Sitter-invariant vacuum is not unique. There are a family of 
solutions corresponding to different vacua [5]. Moreover, the general solutions have 
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two singularities, at z = or Z = 1 (when one point is on the lightcone of the other), 
and at z = 1 or Z = - 1 (when one point is on the lightcone of the point antipodal 
to the other). Among these solutions, however, there is one two-point function that 
is the Wightman function for a unique Euclidean vacuum, or Bunch-Davies vacuum, 
which is only singular when the two points are lightlike related, at z = 0, and which 
can be obtained by analytic continuation from the Euclidean de Sitter space (i.e., a 
four-sphere) , namely [2H [221 ISSl [2H EH [261 El El EZl [28] 

G w (z) = c m 2 F 1 (h + ,h^,2;l - z), (11) 

where h± = (3/2) [1 ± \J\ — (4/9)m 2 ] and 2 F\ is the Gaussian hypergeometric func- 
tion. The constant c m = r(/i + )r(/i_)/(47r) 2 is obtained by requiring the Wightman 
function to have the same singularity strength as that in Minkowski space. Note 
that for x, x' timelike separated, z < 0, one should take care of the branch cut of 
the hypergeometric function by using the ie prescription (c.f. [T8l 129] ) 

, f z(x, x') — ie if t > t' , 1fVl 

Z(X,X') < ; A , • -r . • (12) 

v ' ' I z(x,x') +ie if t' > t K ' 

In the following, the same prescription is also implicit in the logarithmic function, 
where, for timelike separations, 

mz[x,x ) - | ln | z | +i?r if f>t - [id) 

In the Bunch-Davies state, as m — >■ 0, G w in Eq. f JTT]) can be expanded as [6| [T9] 
(with a ^-independent O(m ) term dropped) 

* JL + _L(i_2hi^ +0(m 2 ). (14) 



167T 2 m? 167T 2 

As one can see, the Wightman function is divergent in the massless limit. This 
divergence is eliminated for n-point combinations of the Wightman function that 
are shift invariant. For such shift-invariant operators of a massless scalar field, we 
can drop the mass dependence of the massive Wightman function in the massless 
limit to get 



G = T7~2 ( 7- 21112 ) > ( 15 ) 



16n 2 \z 

the shift-invariant part of the two-point function of a minimally coupled massless 
scalar field in the de Sitter-invariant vacuum state |0), where z = z(x,x') is given 
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by Eq. (|9]) as one-fourth the invariant interval between the two points in the 5D 
Minkowski spacetime in which the 4D de Sitter may be embedded as a unit hyper- 
boloid. 

An equivalent way to get G(x,x') is to start with the spectral representation 
of the two point function Ge{x,x') for a massive scalar field on the Euclidean de 
Sitter space (four-sphere) |30j E] • Note that the Lorentzian Wightman function of 
the Euclidean vacuum can be obtained via continuing the Euclidean z(x,x') on the 
5* 4 to its Lorentzian version, and assuming the same ie prescription as in f fl2|) . Ge 
obeys {D x — m 2 )G j g;(x, x') = —5(x, x') with D x being the Laplacian with respect to 
the x coordinates. Hence this Green's function may be written as 



where the 4> n (x) are an orthonormal set of real eigenfunctions, obeying the equation 
(O x — m 2 )(j) n (x) = —\ n (f) n (x) with eigenvalues X n = m 2 + l n {l n + 3) and having 
(/ + 1)(/ + 3/2) (Z + 2)/3 orthonormal eigenfunctions sharing the same nonnegative 
integer value of l n = I. The lowest eigenvalue, which we shall label as n = with 
Z = 0, is A = m 2 , corresponding to the constant eigenfunction <po(x) = l/y/V^, 
where V 4 = 8vr 2 /3 in our units with H = 1, so that V4 is the volume of the unit S . 

Clearly the n = term in this spectral representation of Ge{x,x') diverges 
when m = 0. One can also directly see that when m = 0, there is no solution 
to (O x — m 2 )G^(x, x') = —5(x,x') on a compact Euclidean manifold such as the 
four-sphere, since then the integral over x of the left hand side is identically zero, 
whereas the integral of minus the covariant delta function on the right hand side 
gives —1. However, if we omit the zero-eigenvalue term in the sum, we get a result 
that is finite even when m = and then is uniquely defined as [301 E] 



Our Lorentzian G(x, x') can thus be regarded as an analytic continuation of G l 
This shift-invariant part of the two-point function obeys the equation 




(16) 



n 




(17) 



n x G%(x,x') = n x ,Gl{x } x') = -5(x,x') + 



(18) 
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so it is not a Green's function for the Laplacian on the four-sphere; such a Green's 
function does not exist. A similar equation holds for the Lorentzian Wightman two- 
point function G(x, x') with the 5 function absent. Therefore, the de Sitter-invariant 
G(x,x') is not the two-point function (if}\<f>(x)^>(x')\if}) in any quantum state of 
the original Fock space, consistent with Allen's proof [5] that there exists no de 
Sitter-invariant Fock vacuum state, in which a de Sitter-invariant two-point function 
would be defined. 

However, in our alternative set of quantum states in which only shift- invariant 
operators are assigned expectation values, the failure of G(x, x') to be a solution to 
the homogeneous equation of motion cancels out. The shift-invariant expectation 
values will never have a single G(x, x') term with the argument x, but terms with 
argument x will always occur in pairs, such as G(x, x') —G(x, x"), and the Laplacian 
with respect to x acting on such a combination will always be zero in the Lorentzian 
spacetime. 



3 Examples of Shift-Invariant Correlation Func- 



We regard the shift invariance as a physical constraint on constructing observables. 
Then there is a de Sitter-invariant vacuum for a massless minimally coupled scalar 
field as long as one considers only shift-invariant operators. In this section, we 
give some examples of these operators and show that their correlation functions are 
indeed both de Sitter invariant and free of IR divergences. 

3.1 Two-Point Functions of Derivatives 

The shift-invariant correlation function for the derivatives of a massless scalar field 
is readily calculated to be 



For two points separated by proper time T = —ifi(x, x'), so that Z = coshT and 
z={l-Z)/2 = - sinh 2 (T/2), we have 



tions 



(O|V a 0(x)V y 0(x')|O) = V a Vy 



G(x,x') 



(4 + 2z)n a ny + (1 + 2z)g aV 
32tc 2 z 2 



(19) 



n a = ie' 




(20) 
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n b ' 



Figure 1: Two timelike separated points x and x', with derivatives of the field along two 
timelike directions. 



where e a and e b ' are unit tangent vectors of the geodesic at x and x' satisfying 
e a e a = e b ' ey — — 1. See Fig. [TJ Then, using Eq. (fT5|) . we have 

3 3 



(0|e»^V^WV^(xO|0> ^eV'V.WGM = ^ - (r/2) 
This correlation function goes to zero exponentially rapidly as T — > oo. 



• (21) 




Figure 2: Two timelike separated points x and x', with derivatives of the field along two 
spacelike directions. 



If e a and e b ' are two unit spatial vectors orthogonal to n a and n b ' , and e 6 ' is 
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obtained by parallel transporting e a along the geodesic, as shown in Fig. |2j then 

l + 2z 1 - 2sinh 2 (T/2) 



(0\e a V a <t*{x)e b \ b 4{x')\0) 



(22) 



32ti 2 z 2 32vr 2 sinh 4 (T/2) ' 

This also vanishes exponentially as T — > oo, though asymptotically at half the rate 
of that of the previous case. 

For two spacelike separated points that are connected by a spacelike geodesic of 
length ji — L, with e a = —n a and e b ' = n b ' = —g ab 'n a now being spacelike tangent 
vectors to the geodesic, we have 

3 3 



(O|e a V a 0(x)e 6 'V 6 ^(x')|O) 



(23) 



32tt 2 z 2 327r 2 sin 4 (L/2)' 

On the other hand, if e a and e b ' = e a g ab ' are not tangent vectors to the geodesic but 
are orthogonal to n a and to n b ' respectively, we have 

1 + 2,2 l + 2sin 2 (L/2) 



(e"V^V^M)- 32A2 - 32i28inl(i/2) . 
3.2 Correlation Functions of Differences 



(24) 



The difference between field operators at different points is shift invariant. So, if we 
calculate the correlation function (cf. [7J Eq. (45) for a similar construction) 



G 



xy,uv 



m— >0 



= lim [G w (x,u)- G w (y, u)-G w (x, v) + G w (y, v)} , 
m— >0 

then the divergent terms as m — > cancel out, leaving only 



(25) 



G 



xy,uv 



G(x, u) — G(y, u) — G(x, v) + G(y, v) 
2 In z(x, u) 

2hxz(x, v) + 



16vr 2 



z(x, u) 

1 



z(y,u) 
1 



— 2 In z(y, u) 

— 2 \nz(y, v) 



, (26) 



z(x, v) 

which is free of IR divergences. This will be further confirmed in the examples given 
later. 

Note that although D x , the Laplacian with respect to x, acting on G(x,x') 
has a nonzero constant term, the combination that appears in G XV)UV does give 
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O x G xy>uv = 0, with the constant term cancelling out. However, the cancellation 
also implies that one could replace the de Sitter-invariant G(x,x') with a suitable 
non-de Sitter-invariant G(x,x') and still obtain de Sitter-invariant shift-invariant 
expectation values such as G XV}UV . For example, G(x, x') = G(x,x') + f(x) + f'(x') 
with arbitrary functions / and /' will give the same G(x, u) — G(y,u) — G(x,v) + 
G(y,v) = G xy . In particular, if one chooses / and /' to obey O x f(x) = 

D x 'f'(x') = — I/V4, then instead of O x G(x,x') = n x iG(x,x') = I/V4, one gets 
O x G(x,x') = n x >G(x,x') = 0, so that G(x,x') obeys the equations of motion of 
a true Wightman function. One might think (as we did) that G(x,x') would be 
the Wightman function of a suitable non-de Sitter-invariant Fock state, but Albert 
Roura [31] has convinced us that this is apparently not the case. 

For example, following Allen [5], one can solve the equation of motion of a 
massless scalar field in the k = FLRW coordinate system of de Sitter spacetime 
(flat spatial slices), 

ds 2 = -dt 2 + e 2t dx 2 . (27) 

Then it seemed that one could choose a vacuum defining all n-point functions (not 
just the shift-invariant ones) that is not fully de Sitter- invariant but only E(3) 
invariant, corresponding to the symmetries of rotations and translations on the flat 
spatial slices. There it appeared that one could choose the two-point function to be 
[5] (up to an arbitrary additive constant) 

D(x,x') = — f- -2\nz + 2t + 2t') . (28) 
lo7r \z J 

Note that this is not de Sitter-invariant due to the last two terms depending on 
the time coordinates, which correspond to f(x) = t/(3Vi) and f'{x') = t'/(3Vi). 
However, in the similarly defined shift-invariant correlation function G xysUV , the t- 
and if:'- dependent terms cancel out. The remaining part is exactly the same as Eq. 
( I26|) . As we see, although the two point function D(x,x') is not de Sitter- invariant, 
it is not shift-invariant either and thus may be regarded as "unphysical." On the 



1 This may be regarded as an analogue to the so-called 'physical equivalence' in [32] , where the 
graviton two point functions A aba > b ,(x,x') and h aba > b >(x,x') = A abalbl (x,x') + V \ a Qb)a>b'{ x i x ') + 
V (a'Q\ a b\b') { x i x ') are physically equivalent in the sense that both give the same two point function 
of gauge-invariant operators, e.g. the linearized Weyl tensor. 
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other hand, the shift-invariant correlation function G xy ^ uv obtained from it is de 
Sitter- invariant . 

However, Allen and Folacci [6], following work by Ford and Vilenkin [33], pointed 
out that the procedure used to obtain the two-point function D(x,x') involved a 
regulator that dropped zero modes, so it is not the two-point function of a Fock 
state, but rather of an "unrealizable limit of a continuous family of Fock vacuum 
states." The two-point functions of actual Fock states have not only the additive 
functions f(x) and f'(x') but also product terms that do not cancel when one forms 
the shift-invariant four-point function G xyjUV from themj 2 ] 

More generally in the k = FLRW coordinate system, one could have f(x) = 
t / (3V4) + C\e~ 3t + C2 with two arbitrary coefficients c\ and C2, and similarly for f'(x'), 
to give O x G(x,x') = O x 'G(x,x') = but the same G XVtUV = G(x,u) — G(y,u) — 
G(x, v) + G(y, v). However, it appears doubtful that there are Fock states that give 
these two-point functions either. 

In the k — — 1 FLRW coordinate system with hyperbolic spatial slices and 
scale factor a(t) = sinht, one can have f(x) = (In sinht + 1/ sinh 2 t)/(3V 4 ) + 
Ci(lntanh(t/2) + cosht/ sinh 2 t) + C2 and similarly for f'(x') to get an analogous 
two-point function G(x,x') = G(x,x') + f(x) + f'(x'). Similar to the case with the 
k = FLRW coordinate system, one can get an idealized limiting state giving this 
two-point function, but the results of |6J suggest that this is also impossible with a 
actual Fock state. 

In the k — +1 FLRW coordinate system with three-sphere spatial slices and 
scale factor a(t) = cosht, one can have f(x) = (lncosht — 1/ cosh 2 t)/(3V^) + 
Ci(arctan sinh t + sinh t/ cosh 2 t) + C2 to get yet another non-de Sitter-invariant two- 
point function G(x, x') = G(x,x') + f(x) + f'(x') obeying D x G(x,x') = D x /G(x,x r ) = 
but the same de Sitter-invariant G xy>uv = G(x, u) — G(y, u) — G(x, v) + G(y, v) as 
obtained from the de Sitter-invariant G(x,x') that does not obey the equation of 
motion for a Wightman function. However, again it appears to be impossible to get 
this G(x, x') as the Wightman function of an actual Fock state [7J [13 l3~Tj . 

2 Indeed, in [B], under the "unrealizable limit" (a — > 0, c.f. Eq. (4.13) of that paper), the 
zero modes of the 0(4) invariant vacuum diverge, but the resulting two-point function is finite 
and essentially corresponds to our G. In particular, the product term disappear and the time- 
dependent terms are exactly the same as our f(x) and f'(x') given below for the k = +1 case 
(after a coordinate transformation). 
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3.3 Some Examples 

One may consider many configurations. One simple example is four points x, y, u, v 
all in order along a timelike geodesic, with x — y and u — v proper time separations 
t, and x — u and y — v separations T > t. Then it is very easy and straightforward 
to show, using the relation and (JSJ), that 

G 1 f 1 



~ 16vr 2 1 S inh 2 (^) ' sinh 2 (^) sinh 2 (f ) 
+ 21n[sinh 2 (^^)sinh 2 (^^)] - 41n[sinh 2 (|)] 
(cosh 2 T + coshT — cosht — 1) (cosht — 1' 



4tt 2 
- In 



(coshT — l)(coshT — cosht) 2 
cosh T — 1 



(29) 



cosh T — cosh t 

and that for T ^> t, G xy>uv — > 0. In the following, we will consider some less trivial 
examples. 

In the static coordinate system, 

ds 2 = -(1 - r 2 )dt 2 + (1 - r 2 Y l dr 2 + r 2 dQ 2 , (30) 

the r-coordinate lines (e.g., lines in which all the other coordinates are fixed) are 
geodesies, i.e. 



V b V h V a = 0, 



where V a = (d r ) a / \[g^r is the normalized unit vector in the r-direction, the deriva- 
tive with respect to the radial proper distance 



f r dr' 

/ = sin -1 r (so r = sin//). (31) 

Jo vl-r' 2 



In general the t-coordinate lines are not geodesic, but an exception is the t- 
coordinate line that passes through the origin. 

The embedding mapping between the 5D bulk space X a (a = 0, . . . , 4) and the 
static patch x^ = (t, r, 6, <fi) is 



X° = VI -r 2 sinh(t), 
13 



X 1 = r sin 9 cos 0, 
X 2 = r sin 9 sin 0, 

X 3 = T COS 



X 4 = Vi -r 2 cosh(t), (32) 

with r & [0,1], 9 £ [0, 7r] and e [0, 27r). Inserting these transformations in Eq. 
gives 

Z(x,x') = — r 2 )(l — r' 2 ) cosh (i — £') + rr' cos fi, (33) 

where 

cos f2 = cos cos 9' + sin sin 9' cos (0 — 0'). (34) 
One can then get z(x,x') from Eq. ([9]). 

3.3.1 Configuration 1 



y 



L 



Figure 3: Configuration 1. (x,y) and (u,v) are regarded as two clusters with a timelike 
separation T. 



For this configuration in Fig. [31 the coordinates of the four points are 

x : {t = 0, r = r L ,9 = 0,0 = 0}, y : {t = 0, r = r L , 9 = vr, = 0}, 

u : {t = t T ,r = r L ,fl = O,0 = O}, v : {t = t T , r = r L , 9 = it, = 0}, (35) 
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where the geodesic distance between x and y (and between u and v) is L, while that 
between x and u (and between y and v) is T. can be determined in terms of T 
via 

z(x, u) = - sinh 2 (T/2) = -(1 - r£) sinh 2 (t T /2). (36) 
The interval between x and y (and between u and v) is 

z(x, y) = z(u, v) = sin 2 (L/2) = r\. (37) 

The interval between x and v (and between y and w) is, using Eq. (1361) and Eq. (l3Tj) . 

zfou) = = i[l + r 2 - (1 - r 2 ) cosh(t T )] = sin 2 (L/2) - sinh 2 (T/2). (38) 



Then the correlation function is 
G 



87T 2 



1 1 / sin 2 L/2 

o ^ 9 + 2 In 1 9 V 7 ; - ie 

sinh 2 (T/2) -sin 2 (L/2) sinh 2 (T/2) V sinh 2 (T/2) 

(39) 

where the ie picks out the branch of the logarithm (with a branch cut along the 
negative real axis) with the imaginary term — 7ti when the argument of the logarithm 
is negative, which is when sinh(T/2) < sin(L/2), so that the two separations of 
(x, v) and (y, v) are spacelike. The shift-invariant correlation function goes to zero 
for T>i, i.e. as the two clusters are separated far away in time. 

3.3.2 Configuration 2 

In this case we choose four points located at 

x : {t = 0, r = r L ,9 = 0,0 = 0}, y : {t = t T , r = r L , 9 = 0, = 0}, 

u : {t = 0,r = r L ,9 = vr, = 0}, v : {t = t T , r = r L , 6 = it, = 0}, (40) 

where the geodesic distance between x and u (and between y and v ) is L, while 
that between x and y (and between v and u) is T, as shown in Fig. HI The relation 
between tx and T can be expressed explicitly by calculating z(x,y): 

z (x,y) = -sinh 2 (T/2) = -(1 - r\) sinh 2 (t T /2). (41) 
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Figure 4: Configuration 2. (x, y) and (u, v) are regarded as two clusters with a spacelike 
separation L. 



The geodesic distance L between x and u (and between y and v) can be charac- 
terized by 

z(x, u) = z(y,v) = sin 2 (L/2). 
Moreover, the interval between x and v (and between y and w) is given by 



z(x,v) = z(y,u) = sin 2 (L/2) - sinh 2 (T/2). 



(42) 



Note that when the ie prescription is explicit, due to the opposite time ordering of 
(x, v) and (y, u), z(x, v) and z(y, u) are not equal but differ by the sign of ie. There- 
fore, even when these two separations are timelike, their imaginary contributions to 
the Wightman functions cancel. 
Then the correlation function is 

sinh 2 (T/2) 



G 



xy,uv 



1 



2 In 



1 



sin 2 (L/2) 



(43) 



_sin 2 (L/2) sin 2 (L/2) - sinh 2 (T/2) 

which goes to zero as T L. The absolute value here is due to the proper treatment 
of the ie, with one +ni canceling another —iri when sinh(T/2) > sin(L/2). Kirsten 
and Garriga [7] have previously noted the linear divergence in the time as T — > oo. 
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4 Discussion 



We have shown that if one restricts the definition of a quantum state to giving 
expectation values of shift-invariant quantities (such as products of differences of 
field values) involving a massless scalar field, then there does exist a perfectly well- 
behaved de Sitter- invariant vacuum state for a massless scalar field. It just does not 
give well-defined expectation values to quantities that are not shift invariant, such 
as the product of two field values. Since the Lagrangian of a massless scalar field is 
shift invariant, it is natural to restrict the expectation values given by the quantum 
state to shift-invariant quantities, analogous to the way that one restricts to the 
expectation values of gauge-invariant quantities for a gauge-invariant Lagrangian. 

Our results were perhaps first partially anticipated by Pathinayake, Vilenkin, 
and Allen [34] , who showed that the related theory of a massless antisymmetric 
tensor field (a two-form potential with the Lagrangian given by the trace of the 
square of the three-form field strength that is the exterior derivative of the two- 
form) does have a de Sitter-invariant quantum state with a well-defined two-point 
function. When the equations of motion are satisfied, the dual of the field strength is 
the gradient of a massless scalar field that obeys its equations of motion. Therefore, 
the expectation values of combinations of the dual of the three-form field strength 
can be interpreted as giving the expectation values of gradients of a massless scalar 
field. In principle one can take combinations of line integrals of the dual of the field 
strength to give differences of the scalar field at the two endpoints of the lines and 
hence to get the expectation values of the shift-invariant operators of the massless 
scalar field considered in our paper, but that was not done explicitly in [34J. 

Folacci [30] gave similar results as ours by a BRST quantization on the Euclidean 
four-sphere S 4 but did not consider restricting the quantum state to giving only the 
expectation values of shift-invariant operators on the Lorentzian de Sitter spacetime. 
He therefore concluded that the infrared divergence in Lorentzian de Sitter is real. 

Bros, Epstein, and Moschella [9j, following earlier work in two-dimensions by 
Bertola, Corbetta, and Moschella pj, also used a construction similar to ours of 
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removing the Euclidean zero-eigenvalue eigenfunction (the constant function) and 
then noting that this gives "a local de Sitter invariant quantization of that field on 
the space of test functions having zero mean value." As Moschella expresses it |35J, 
this present paper "is exploring concrete states belonging to the physical subspace 
we have generally constructed." 

Kirsten and Garriga [7] have given the results perhaps most obviously similar 
to ours, a de Sitter-invariant state that is not normalizable. They note, "This 
should not be taken as an indication that the state is pathological: it simply means 
that all values of [the spatial mean of the scalar field] are equally probable." A 
Euclidean derivation of this de Sitter-invariant state has been given by Tolley and 
Turok [36]. What is perhaps somewhat new in our present paper, besides the detailed 
results, is the explicit description of the quantum state as being defined by giving 
the expectation values of just the shift-invariant operators. 

In our detailed results we have shown, in the case of field derivatives, that the 
two-point function becomes vanishing as the timelike separation goes to infinity. In 
the case of field differences, the two clusters (x, y) and (it, v) in both Configuration 
1 and 2 become uncorrelated as the separation between the two pairs of points 
becomes large, although there is a divergence in Configuration 2 in which one takes 
the product of differences of fields between points that are moved apart to arbitrarily 
great separation. 

In other words, the shift-invariant expectation values obey a cluster-decomposition 
property in that the field derivatives become only weakly correlated at widely sep- 
arated points. For the product of field differences, if each of the clusters in which 
the field differences are taken is kept at fixed size but separated widely from the 
other cluster whose field differences are multiplied by the first difference, then the 
expectation value goes to zero. That is, the correlation between the differences of 
fields at two point-pairs that each have fixed separation goes to zero as one pair is 
widely separated from the other pair. There is no IR divergence for such quantities, 
but only for products of differences of fields in which the separations of the points 
whose different field values are used in the product are taken to infinity. 



18 



Although the massless scalar field is only a poor-man's toy model for gravity, it 
does suggest that there is a de Sitter-invariant gravitational quantum state with no 
IR divergences for the correlations of local gauge-covariant quantities, such as corre- 
lations between the Weyl tensor at one location and at another. Mora and Woodard 
[15], and Mora, Tsamis, and Woodard [IB] , correcting some minor errors in earlier 
work by Kouris [2], have recently confirmed this at the linearized level, though, 
contrary to our own opinion, they conjecture that this lack of IR divergences would 
not be true at the next order. Cf. similar results by Frob, Roura, and Verdaguer 
that include the effect of one-loop corrections from matter fields [37] and apply for 
the full Riemann tensor [35] . 

On the other hand, the divergence one gets for the expected product of differ- 
ences of massless scalar field values when the points giving the differences are pulled 
apart to infinity may be analogous to the divergences in the gauge-fixed graviton 
propagator when the gauge fixing is taken over a spacetime separation that is taken 
to infinity (see, e.g., [101 |39j HQ]). For example, in a classical model for de Sitter 
spacetime with quantum fluctuations that is a classical spacetime with small metric 
perturbations, one may define the following two-time function along a central time- 
like geodesic that might well diverge in the limit that the proper time between the 
two times goes to infinity: 




Figure 5: Distortion of a small 'sphere' evolving along a central timelike geodesic in 
perturbed de Sitter spacetime. 
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As illustrated by Fig. [5j if t is the proper time along the central timelike geodesic, 
let the two times be t = and t = T. At the point t — along the central timelike 
geodesic, construct the full set of spatial geodesies orthogonal to the timelike geodesic 
and form the locus of points at proper distance r from the t = point, a two- 
dimensional 'sphere' of proper radius r. Parallel-propagate the four-velocity of the 
central timelike geodesic along each of the spatial geodesies to the 'sphere' and 
use these timelike vectors as the four-velocities of new timelike geodesies running 
forward in time from the 'sphere' at t = 0. Go forward by a proper time T along 
each of these new timelike geodesies to get a second approximate sphere that would 
be nearly simultaneous to the t = T point along the central timelike geodesic. 

Now construct the following dimensionless invariant measures of the distortion 
of the shape of this second approximate sphere: a bulk 'ellipticity' E that is the 
variance of the proper radial distance from the t = T point to the sphere, divided 
by the square of the average radial distance, and a surface 'distortion factor' F 
that is the area of the approximate sphere, multiplied by the integral over the area 
of the square of its intrinsic Ricci scalar, minus the value this product has for a 
round sphere, which is 64ir 2 . For pure de Sitter spacetime, the linear size of the 
initial sphere (if it is small, r < 1) will just grow by a factor of coshT, but E 
and F will remain zero. For distortions of the spacetime with wavelength much 
longer than the size of the sphere, and for a small sphere size, r < 1, the initial 
'sphere' will distort to an approximate ellipsoid, and for small distortions, E and 
F will be approximately quadratic in the normalized differences in the semimajor 
axes of the approximate ellipsoid and hence approximately quadratic in the metric 
perturbations in a synchronous gauge. 

Given the t = and t = T points that can be taken to define the central 
timelike geodesies (at least in the absence of caustics), and given the initial radius 
r, the measures E and F of the distortion of the approximate sphere at time t = T 
obtained by evolving the initial 'sphere' of proper radius r at time t = are perfectly 
gauge-invariant quantities, but since the timelike geodesies of proper time length T 
that go into the definitions can be continually distorted by the fluctuating spacetime, 
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these T-dependent gauge-invariant quantities E and F need not remain finite as T 
is taken to infinity. They are not localized quantities, or products of two localized 
quantities (like the product of two Weyl tensors at different points), so even in a de 
Sitter-invariant quantum state in which E and F do not depend on the initial t — 
point or the tangent vector there of the central timelike geodesic, these nonlocal 
gauge-invariant quantities E(T) and F(T) can in principle increase without limit 
as T is taken to infinity. This does not mean that the de Sitter-invariant quantum 
state has local fluctuations that are in any sense growing from one location or time 
to another (which would violate de Sitter invariance) but just that gauge-invariant 
quantities like E(T) and F(T) that probe a whole region of spacetime, which gets 
larger and larger the larger one makes T, can grow indefinitely with T, analogously 
to the way that the product of differences of massless scalar field values can grow 
with the proper time separations between the points whose difference of field values 
are used in the product. 

In conclusion, we have shown that there is a de Sitter-invariant quantum state 
for a massless scalar field in the de Sitter spacetime, so long as one just requires the 
state to give well-defined expectation values for operators that are shift-invariant un- 
der constant shifts of all the scalar fields used in the operator. This is a poor-man's 
analogue for a de Sitter-invariant quantum state of gravitational field fluctuations 
that are just required to give well-defined expectation values for gauge- invariant 
operators. If these operators are products of localized gauge- covariant operators in 
two regions (like the Weyl tensor at two points), the cluster-decomposition property 
exhibited for the massless scalar field suggests that the expectation values of the cor- 
responding gauge-invariant gravitational field operators will also tend to constants 
(presumably zero for quantities that vanish in classical de Sitter spacetime, such as 
products of Weyl tensors) in the limit that the two regions are pulled far apart, but 
if the operators involve a whole spacetime region of time period T, it would not be 
surprising for them to diverge in the limit that T is taken to infinity, just as the 
product of two scalar field differences does when the separations of the points whose 
fields are differenced is taken to infinity. 
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